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The classical Hall effect presents a surprisingly unusual and challenging problem in electrostatics, 
with boundary conditions that are not of Dirichlet, Neumann, or of mixed Dirichlet and Neumann 
type. These unusual boundary conditions create several difficulties not normally encountered in 
standard problems, and ultimately lead to expansion of the electric potential in a nonorthogonal 
basis set. We derive the boundary conditions for the potential in a rectangular geometry, construct 
a solution for the potential, and discuss the relation between this problem and problems of the 
standard mixed type. We also address a commonly encountered misconception about the current 
distribution. © 1998 American Association of Physics Teachers. 


I. INTRODUCTION 

The Hall effect was discovered over 100 years ago, and 
has since become a widely used experimental tool for study- 
ing the transport properties of materials, as well as the basis 
for a large number of technological applications. 1 One does 
not need to know the full electrostatic solution to the Hall 
problem in order to extract useful information, since voltage 
differences between suitable pairs of points in the current 
flow suffice to characterize the transport processes. 2 The full 
solution to the electrostatic problem is a surprisingly chal- 
lenging exercise, going beyond the techniques that are most 
frequently used in potential theory problems. 

A typical arrangement of a Hall effect experiment is illus- 
trated in Fig. 1. A thin, rectangular metal plate lies in the 
x—y plane. The edge at y = H is maintained at electric po- 
tential V Q and the edge at y=—H is maintained at — V 0 . 
Thus an electric current will flow in the — y direction. If an 
externally produced uniform magnetic field B is imposed in 
the — z direction, there will be a magnetic force on the mov- 
ing (positive) charges, directed in the +x direction. This 
gives rise to a charge separation that produces a potential 
gradient in the x direction. The problem, then, is to find an 
analytical expression for the electric potential V(x,y) every- 
where in the metal plate. 

Although this problem is stated in simple terms, it turns 
out to have several unusual features. First, it is not immedi- 
ately clear what conditions must be imposed on V(x,y) at 
the boundaries x = 0 and x = L. Second, once these boundary 
conditions are derived, they result in a boundary-value prob- 
lem that is not of standard type. Usually we expect to en- 
counter a Dirichlet problem (V specified everywhere on the 
boundary), a Neumann problem (normal derivative of V 
specified everywhere on the boundary), or a “mixed” prob- 
lem (Dirichlet conditions on part of the boundary and Neu- 
mann on the remainder). 3 As we shall see, our simple prob- 
lem falls into none of these classes. Third, as a consequence 
of these unusual boundary conditions we are forced to ex- 
pand V{x,y) in terms of nonorthogonal basis functions. Nev- 
ertheless, a series solution for V(x,y) can be obtained by 
straightforward methods. This combination of unusual math- 
ematical features and a simple physical situation will, we 
hope, make the problem of some interest to readers of this 


Journal. Finally, as we show in an Appendix, this problem 
can also be solved by reducing it to an infinite number of 
problems of standard type. 


II. BOUNDARY CONDITIONS AT x = 0 AND x = L 

We would like to determine the electric potential V(x,y) 
everywhere in the conductor. In the steady state there is no 
volume charge density so V -E = 0 (E is the electric field). 
As this is electrostatics, E=-Vk, so we need to solve 
Laplace’s equation, 

v 2 y=o, (l) 

subject to the appropriate boundary conditions. The upper 

and lower edges of the plate are maintained at constant po- 
tential: 

V(x,H) = V 0 , (2) 

V(x,-H) = -V 0 . (3) 

The other physical constraint is that no charge enters or 
leaves the sample through the sides at x = 0 and x = L. We 
must now express this condition on the electric current in 
terms of the electric potential to provide us with the bound- 
ary conditions for the left- and right-hand edges. 

In the case of zero magnetic field, the current density J 
and the electric field E are related by Ohm’s law (throughout 
this analysis we assume linear materials) 

E=pJ, (4) 

where p is the resistivity of the material. Equation (4) repre- 

sents a condition of balanced forces. The left side is the 
electrostatic force per unit charge and the right side is the 
negative of the drag force per unit charge. In the presence of 
a magnetic field, B, the magnetic force per unit charge, vXB, 
must be added to the left side of Eq. (4). The current density 
is related to the drift velocity, v, as J = nq\, where q is the 
charge and n is the number density of the charge carriers, 
assumed to be constant. Thus the generalization of Eq. (4), 
valid with a magnetic field present, is 

1 

EH JxB = pJ. (5) 

nq 
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Fig. 1 . Sample geometry showing the axes, top and bottom boundary con- at t * le J i s straight down since J x - 0. (c) at the top edge E is 

ditions, and the magnetic field into the page. The points labeled P and Q are straight down since v — H is an equipotential. 

symmetry points. 


Let us consider the standard problem, in which the mag- 
netic field is perpendicular to the plane of the conducting 
sheet: B= — B Q z. In general, within the sample, J may have 
both x and y components. Evaluating the cross product in 
Eq. (5) we obtain 


B o 

E. — pj , J = 0, 

1 nq y 

(6) 

Bq 

E v -pJ v + — J = 0. 
J y nq 

(7) 

Solving Eqs. (6) and (7) for the components of the current 
density we have 4 

(■ E x ~\E y ) 

x p( 1+xV 

(8) 

(E y + \E X ) 
y P( 1+X 2 )’ 

(9) 

where X = B 0 / pnq. 

The requirement that no current leave through the edges at 
x = 0 and x = L means that J x (0,y) = 0 and J x (L,y) = 0. 
These conditions, along with Eq. (8), result in 

E x (0,y) = \E y (0,y), 

(10) 

E x (L,y) = XE y (L,y). 

(ID 

In terms of the electric potential we have E x = —dVIdx and 
E v = — dV/dy, so at both x = 0 and x = L we require 

dV dV 

dx ^ dy ’ 

(12) 


This is the desired boundary condition on V(x,y) at the left 
and right boundaries. Note that the slopes of all the equipo- 
tential curves are equal to — X at the places where they meet 
the left- and right-hand edges. 


The parameter X is also related in a simple way to a pa- 
rameter called the Hall angle 5 d H , defined as the angle be- 
tween J and E. From the vector triangle 6 of Eq. (5), illus- 
trated in Fig. 2, we see that the Hall angle is simply 

0 H =tan~ 1 (B o l pnq) = tarT l X. (13) 

The magnitudes and directions of J and E may well be dif- 
ferent at different places in the plate. Nevertheless, the angle 
0 H between J and E is everywhere the same. 

The vector diagram provides a more physical way of un- 
derstanding the boundary conditions. Figure 2 illustrates the 
vector triangle at (a) an arbitrary position in the rectangle, (b) 
a point on the right edge, and (c) a point on the top edge. On 
either the left- or right-hand edges J must be parallel to the 
edge, thus E makes an angle 0 H with respect to the boundary 
[Fig. 2(b)]. This is precisely the condition expressed by Eq. 
(11). The top edge is an equipotential, so E must be perpen- 
dicular to the boundary [Fig. 2(c)], 

Note that a statement frequently encountered in textbooks 7 
stands in need of a correction. It is often asserted that in the 
steady state the magnetic force on the charge carriers just 
balances the horizontal component of the electric force 
(which is due to the Hall potential) and that, consequently, 
the current flows parallel to the y axis. This is so only at the 
left and right edges (as in Fig. 2). Elsewhere in the plate, J in 
general has an x component and the magnetic force is not 
entirely in the x direction. However, for the case of small 
magnetic field or long, thin samples, it is a reasonable ap- 
proximation. We shall investigate the current distribution in 
Sec. VI. 

III. STATEMENT OF THE BOUNDARY- VALUE 
PROBLEM 

The mathematical problem is as follows: Find V(x,y) that 
solves V 2 V=0 in the region O^x^L and — H^y^H sub- 
ject to the boundary conditions: 

V(x,H) = V 0 , (14) 
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V(x-H)=-V 0 , 


( 20 ) 


dV 

dx 


x=0 


dV 
=X —I 

dy 


r = 0 


dV 

dx 


= X 


x = L 


dV 

dy 


x = L 


(15) 

(16) 

(17) 


Note that A. depends on B 0 so the “no magnetic field” case 
corresponds to \ = 0. 

This is an unusual set of boundary conditions. Typically, 
either the potential (Dirichlet), or the normal derivative of 
the potential (Neumann), is known on each boundary, or else 
the problem is of the “mixed” type, with Dirichlet condi- 
tions on part of the boundary and Neumann on the remain- 
der. In our case the boundary conditions for the left and right 
edges are given in terms of both partial derivatives of the 
unknown potential, V{x,y). Note that the problem reduces to 
the standard “mixed” type for \ = 0. The general (\A0) 
problem does not appear to be covered in the standard treat- 
ments of potential theory. 8 


IV. SOLUTION 

We seek solutions of Laplace’s equation that satisfy con- 
ditions (16) and (17). Any linear combination of such solu- 
tions will still satisfy (16) and (17). The appropriate linear 
combination can then be built up to satisfy Eqs. (14) and 
(15). We will handle the two forms for solutions of Laplace’s 
equation separately. 


d 2 Y 
dy 2 


= k 2 Y , 


where the separation constant, k 2 , is yet to be determined. 
X(x ) and Y (y) are, of course, real functions of their argu- 
ments. 

We can easily show that k must be real, and therefore X is 
a trigonometric function of x and Y is an exponential func- 
tion of y. To see this put V=XY in Eq. (16) to find 


X' 0 Y=XY’X 0 , (21) 

where the prime denotes differentiation with respect to the 
argument (X' =dX/dx,Y' =dY/dy). The subscript 0 de- 
notes evaluation at x = 0. Differentiating Eq. (21) with re- 
spect to y we get 


Xq Y' = XY"X 0 . (22) 

Substituting Eq. (20) in Eq. (22) we find 

X' (l Y' = \k 2 YX Q . (23) 


Then using Eq. (21) in Eq. (23) to eliminate Y' we have 



(24) 


Consequently, k 2 is real and non-negative. Therefore, k is 
real and from the forms of Eqs. (19) and (20) we see that 
X{x) is trigonometric and Y (y ) is exponential. 

For a given value of k, the harmonic-exponential solution 
is then of the form 


A. Linear solution 


V k {x,y) = {A k cos kx + B k sin kx)(C k e ky + D k e ky ). (25) 


A linear function of x or y will clearly be annihilated by 
the Laplacian operator. The bilinear solution V(x,y) = (a 
+ bx)(c + dy), with a,b,c,d constant, satisfies the differen- 
tial equation. However, it is easy to see that this form cannot 
satisfy the boundary conditions Eqs. (16) and (17). Similarly 
V(x,y) = (a + bx)f(y) or (c + dy)g(x), where / and g are 
arbitrary functions, cannot satisfy Eqs. (16) and (17). We are 
therefore left with the possibility of a linear function of x 
plus a linear function of y. At the left or right edge, the slope 
of the equipotential curve is — X. Thus the linear solution 
must be of the form 

V{x ,y) = aX(x + b) + ay + c , 

where a,b,c are constants. We know from elementary con- 
siderations that in the absence of a magnetic field the solu- 
tion is V{x,y) = Voy/H. The case of “no magnetic field" 
corresponds to X = 0 so a = V ( , / 1 1 and c = 0. The linear so- 
lution is therefore of the form 

V(x,y)= [\(x + /?)+y]. (18) 

B. Harmonic-exponential solutions 

For the nonlinear forms we assume separable solutions of 
the form V(x,y) =X(x)Y(y). In the usual way, Laplace’s 
equation then separates into two ordinary differential equa- 
tions: 

d 2 X 


C. Applying the boundary conditions at x = 0, L 

We now apply conditions (16) and (17) to determine the 
separation constant k and to restrict the range of possibilities 
for the constants of integration A, B,C,D. First, applying Eq. 
(16) to Eq. (25), we find 

B k C k e ky -D k e~ ky 

i= x c> +D ;,-v <“> 

Thus the right side of Eq. (26) cannot be a function of y. 


Hence, either 


C k = 0 , B k = Ic ■> 

(27) 

or 


D k = 0, B k = XA k . 

(28) 


For a specific value of k, the most general solution is a linear 
combination of both possibilities. Thus we have 

V k (x,y) = R k {cos kx + X sin kx)e ky 

+ S k ( cos kx — X sin kx)e~ ky , (29) 

where R k =A k C k (for the case D k = 0) and S k =A k D k (for the 
case C k = 0). 

If we impose condition (17) on Eq. (29) and note that the 
terms in e ky and e ky must satisfy this condition separately, 
we find 

(1 +\ 2 )sin kL = 0. (30) 

Hence sin kL = 0 and k is restricted to values k n given by 
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k„ = nv/L, (31) 

where n = 1,2,3,... . 

The most general solution satisfying the boundary condi- 
tions at the left and right edges is therefore Eq. (18) plus 
expressions of the form of Eq. (29), with k„ restricted by Eq. 
(31): 


with the same magnitude but opposite sign. Imposing this 
symmetry on Eq. (33) requires that b = — LI 2. Furthermore, 
both the cosine and sine terms in the sum of Eq. (33) must 
satisfy the symmetry separately. (Remember that X. is a free 
parameter.) Now, since k n = mr/L , 

cos k n (L — x) = (— 1)" cos k n x. 


V(x,y)= [ X (•*' + /?) +y] 

+ 2 [R„(cos k n x+\ sin k n x)e k " y 

n= 1,2,3,... 

+ S„(cos k n x—\ sin k n x)e ~ k " y ] (32) 

or, recombining terms, 

y 

V(x,y)= — [\(x + b) + y\+ 2 [cos k n x{R n e k " y 

n- n= 1,2,3,- 

+ S n e~ k " y ) + \ sin k„x(R n e kny — S„e~ kny )~\. (33) 

D. Exploiting a symmetry 

To proceed further, let us note a symmetry of the system 
of Eqs. (1), (14), (15), (16), and (17). These have the prop- 
erty that V — > — V when x—>L—x and y— + — y. This means, 
for example, that points P and Q in Fig. 1 have potentials 


sin k n (L—x) = {— 1)" + 1 sin k n x. 

Requiring that each term in Eq. (33) satisfy the symmetry 
V— > — V when x—xL — x and y— » — y then leads to (from the 
cosine terms) 

(-1 ) n (R n e- k » y + S n e k » y )=-R n e k » y -S n e- k » y , 

and to (from the sine terms) 

(-1 r + 1 (R„e- k ' y -S„e k "y)=-R n e k *> + S n e- k *y. 

Adding these two equations, we find 

S n = (-l) n+1 R n . (34) 

Two cases arise: 

for odd n , S„ = R„ , 

(35) 

for even n, S n =—R n . 

Using Eq. (35) in Eq. (33) we arrive at our final form of the 
solution: 


V 0 

V(x,y)= — 


X(v--|+y 


+ E T„ 

m — 1,3,... 


cos[ —j— x I coshl —j— y I 4- A. sinl —j— x I sinhl -j— y 


+ 2 u, 

n = 2,4,... 


cos[ —j~ x I sinhl —j— y I + \ sinl —j~ x I coshl —j— y 


(36) 


where T m = 2R m (for m odd) and U n = 2R n (for n even). 
From here on we shall use the subscript m to label the odd 
coefficients T x , T 3 , . . . and the subscript n to label the even 
coefficients f/ 2 ,f/ 4 ,... . 

E. Applying the boundary conditions at y = ±H 


Adding and subtracting the two versions of Eq. (37) gives 

V 0 


H 


L 

Mx~- 


= 2 T n 
m = 1,3,... 


m 77 \ m 7 T 

cosl — — x cosh — — H 


We now require that the T m and U„ be chosen so that 
V (x,y) satisfies the top and bottom boundary conditions 
V(x,±H) = ± Vq . Thus we require 


V 0 

±y ° = 7/ 


X| x--\±H 


+ S T„ 

m = 1,3,... 


±\ sin 


+ 2 u n 

n = 2,4,... 


1 m 77 \ 

rmr \ 

cosl —j— x | cosh 

(- h ) 

\ / til 77 \ 


x | sinhl —j— H \ 


In 77 \ 

HTT \ 

± cos x sinh — — H 

i i 

\ L 1 


I n 77 n it 

+ X sin x cosh H 
\ L L 


(37) 


+ X 2 


u„ 



n 7r 

cosh —— H 
L 


and 


(38) 


rnn mu 

0 = X 2j T m s i n ~T~ x s ' n h| — — H 

m = 1.3.... \ L I \ L 


x , / n 77 \ n 77 \ 

+ 2j U n cosl — x I sinhl — H I . (39) 

In Eqs. (38) and (39) we have on the left-hand side a func- 
tion of x and on the right side its expansion in a series of 
sines and cosines of multiples of ttx/L. It is important to 
note that these functions of x do not form an orthogonal 
basis on the interval O^x^zL. The functions coslnnrx/L) 
and cos(m' itx/L) are orthogonal (for m=f=m') on the interval 
O^Zx^sL. Similarly, sin (mrx/L) and sinOr ttx/L) are orthogo- 
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nal. But sm{mrxlL) and cos/nnrx/L) (with n even and m odd) 
are not. These functions are orthogonal on the interval 
0=£x^2 L, but are not orthogonal on the interval O^x^L. 
This is another unusual feature of this simple Hall effect 
problem. The nonorthogonality of the basis functions results 
from the boundary conditions Eqs. (16) and (17). We must 
be careful in evaluating the expansion coefficients. 

(To understand the nonorthogonality an analogy might be 
helpful. Consider writing a vector in terms of a particular set 
of nonorthogonal basis vectors and then determining the 
components. Let a vector in the plane. A, be written in terms 
of a pair of unit vectors u and v which are at an angle a. We 
can write A =A u u +A v v. To find the coefficient we take the 
dot product of A with the appropriate basis vector. For ex- 
ample to find A u we have u-A =A u u-u+A v ii-v=A u 
+A V cos a. We do not get just the coefficient A u is now 
mixed in as well. The presence of the other coefficient is due 
to the lack of orthogonality. For the usual case with orthogo- 
nal basis vectors a=7r/2 and the second term vanishes.) 

Multiply Eq. (38) by cos(m' ttx/L) and integrate from x 
= 0 to x = L, with the result 

4 V 0 LX 

T — _ 

m 7r 2 Hm 2 cosh(/H7r///L ) 

4X In 7 tH\ n 

tt cosh(nnrH/L) n = 2 A,... ^ n C ° S \ L ) n 2 — m 2 ' 

(40) 

Multiply Eq. (39) by cos(n ' ttx/L) and integrate from x = 0 to 
x = L to obtain 


through order X 2 , which is certainly adequate for comparison 
with most experimental results. 


G. Solution to any order in X 


The solution to any order in X may be obtained by itera- 
tion with Eqs. (40) and (41). A more systematic approach is 
as follows. Let T denote a vector formed from the T m , and U 
a vector formed from the U n : 


/ M 


/ 0 \ 

0 


u 2 \ 

77 


0 

0 

u 4 

\ T :l 


:/ 


Then Eqs. (40) and (41) may be written 
T = T°+ ©U, 


U=17T, 


where the components of T° are 


T 


o 

m 


4 V 0 L\ 

TT 2 Hm 2 cosh(imrH/L ) 
0 otherwise. 


for m odd 


The matrix elements of © are 


(44) 


(45) 

(46) 


(47) 


© 


mn 


_ — 4X „ 

^ n 77 sinh(» ttH/L) m = 1,3,... " 


sinhl 


m ttH 


,2 • 


(41) 


(" 4 X cosh(mrH/L) 

77 cosh(w ttH/L) n 2 


for m odd and n even 


[ 0 otherwise. 


These are our conditions on the coefficients, 7j„ and U„ . 
Note that they involve coupled infinite sums. The odd coef- 
ficients, T m , are expressed in terms of the even ones, U n , 
and vice versa. This results from the failure of our basis 
functions to be orthogonal. 

The solution for the potential V{x,y) is therefore given by 
Eq. (36) with T m and U n determined by Eqs. (40) and (41). 

F. Solution to order X 2 

For typical metals in experimental situations 9 X 
*“0.005-0.1. This suggests that Eqs. (40) and (41) can be 
solved iteratively. Looking at Eq. (40) we see that to first 
order in X we need only the first term. The second term has 
an explicit X and depends on the U n , all of which involve X. 
Therefore, to first order in X the T m are 

4 V 0 LX 

t ~ (421 

m TT 2 Hm 2 cosh(m ttH/L) 

Substituting Eq. (42) into Eq. (41), we find 


The potential V(x,y) is therefore given by Eq. (36) with 7j„ 
and U „ given by Eqs. (42) and (43). This solution, with the 
coefficients cut off as in Eqs. (42) and (43), is correct 


— 16 V q L\ 2 „ Vdx\\\(m ttH/L) 

tt 3 H sinh(n ttH/L) m(m 2 — n 2 ) 


(43) 


(48) 

and the matrix elements of 17 are 


n 


nm 


— 4X sinh(w ttH/L) m 
77 smh(mrHIL) m 2 — n 2 
0 otherwise. 


for m odd and n even 
(49) 


Substituting Eq. (46) into Eq. (45), we have 
T = T° + ©17T 


(50) 


or 

(7-©n)T=T°, (51) 

where I is the identity matrix. Thus the solution for T is 

T=(/-@a) _1 T°, (52) 

T=(/+@a + (©f7) 2 + (©a) 3 + ---)T°. (53) 

Now, as Eqs. (48) and (49) show, both © and 17 are propor- 

tional to X. Equation (47) shows that T° is also proportional 
to X. Hence T only contains elements with odd powers of X. 
And thus, by Eq. (46), U contains only even powers of X. 
The truncated solution (to order X 2 ) given above, with the 
T m and U n given by Eqs. (42) and (43) is equivalent to 

T— T°, (54) 
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Fig. 3. Equipotential curves over the whole conductor. From top to bottom 
the contours shown correspond to V= 1.0,0.95,...,— 0.95,— 1.0. The potential 
was determined as discussed in the text with A. = 0.2, Vo= 1, H=l, L— 1. 


y 



Fig. 4. Equipotential curves for the upper portion of the sample; this is an 
enlargement of Fig. 3 with the same parameters. Top to bottom the contours 
shown correspond to V= 0.99,0.98,0.97,... . 


no magnetic field the equipotentials are horizontal lines. 
With a magnetic field the equipotentials are clearly curved. 
The equipotentials get closer together as one approaches the 
upper left and lower right corners. This corresponds to a 
stronger electric field in these regions. However, as expected 
from the boundary conditions, the equipotentials intersect the 
left and right edges at a constant angle, independent of y. 
This can be seen in the close-up of a portion of the conduc- 
tor, corresponding to the upper part, 0.8^y^ 1, as shown in 
Fig. 4. From top to bottom the contours correspond to V 
=0.99, 0.98, 0.97 

VI. THE CURRENT DISTRIBUTION 


U=UT°. (55) 

Although for real experimental situations A. is a small param- 
eter, in principle one may solve the problem by matrix in- 
version, using Eq. (52), even for large In this case, one 
must truncate the matrices to a finite number of elements in 
advance. 


With Eq. (36) in hand for the potential, it is straightfor- 
ward to investigate the current distribution. Using Eqs. (8) 
and (9) and taking appropriate partial derivatives of V to get 
the components of E we find 

1 hi 77 

J x = — 2 j T m —— sm(imrx/L)cosh(nnry/L) 

P m= 1,3,... ‘ 


V. NUMERICAL RESULTS 

To get a physical sense for our solution we evaluate the 
potential and the current density numerically. The solution is 
given by Eq. (36) and the expressions for the coefficients, T n 
and U m , are given by Eqs. (52) and (46), respectively. As 
the coefficients are infinite in number we truncated the vec- 
tors by taking 20 elements in both T and U and the corre- 
sponding © and <1 were 20X20 matrices. The sums in Eq. 
(36) went to m = 19 and n = 20, respectively. Reasonable val- 
ues of the parameters were chosen to show representative 
behavior. In what follows we let A. = 0.2, V 0 = 1, H= 1, and 
L = 1. These last three parameters just set the scale and are in 
arbitrary units. 

Equipotential curves, V= const, for the conductor are 
shown in Fig. 3 for V= 1,0.95,..., — 0.95,— 1. For the case of 


1 n tt 

H — 2j U n —— sm(mrx/L)sinh(mry/L), (56) 

P n = 2,4.... T 


y Hp 


1 


j. ill. n 

2j T m — cos(OT7rx/L)sinh(;HT7y/L) 

P m = 1,3,... L 

1 n tt 

2j U n —r~ cos(«77x/L)cosh(n'7ry/L). 

P « = 2,4,... L 


(57) 


It is possible to perform two easy tests to verify that J as 
calculated is consistent with the original boundary condi- 
tions. 
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First, and most obviously, J r vanishes at x = 0 and x = L. 
This is consistent with the boundary condition illustrated by 
Fig. 2(b). 

Second, we shall show that at y=±H, J x = —\J y . This 
means that at the top or bottom edge, the current makes an 
angle 0 H = tan -1 X with the normal to the boundary, as 
shown in Fig. 2(c). That this condition is met by Eqs. (56) 
and (57) is not immediately obvious. Therefore, differentiate 
Eq. (37) with respect to x and divide by p to obtain 


Va\ x 1 m 77 

°=-^+ 2 - T m — 

pH m = 1 . 3 .... P L 


— sin 


x cosh 


H 


m 7r 

± X cos[ —j— x | sinh 


niTT 


H 


x , 1 n 77 

+ 2 - U n ~j~ 

n = 2 , 4 ,... P ^ 


h- sin — — x sinh — — H 


n 77 \ / n 77 

+ X cos| —j— x | cosh( — — H 


(58) 


Evaluate Eq. (56) at y=±H , and to the resulting equation 
add Eq. (58) to obtain 


VaX X mir 

J x (x,±H )= — — ± — 2j T m — — cos(rmrx/L) 

pH p tn = 1,3,... f 


X n tt 

X sinh( m ttHIL ) H — 2j U n —r~ 

P (1 = 2 , 4 .... L 

Xcos(mTx/L)cosh(mrH/L). (59) 

Comparing Eq. (59) with Eq. (57) evaluated at y = ±H, we 
see indeed that 

J x (x,±H)= — \J v (x,±H) (60) 

for all x. Thus all across the top boundary, the current enters 
at an angle Oh from the normal. 

It is worthwhile to note one final point about the y com- 
ponent of the current. Integrating Eq. (57) across a horizontal 
section gives the net current in the y direction: 

[L ~ V Q Lt 

Iy=t )o JydX= ~pH~’ (61) 


where t is the thickness of the conducting plate. This is just 
the current we would expect in the absence of a magnetic 
field. This is why the standard elementary treatments, which 
ignore the existence of a transverse current, nevertheless ob- 
tain reasonable behavior. Although the current is not every- 
where in the y direction, the net y current does not depend 
on the magnetic field. 

To obtain numerical results for the current density we use 
the same procedure for determining the coefficients T m and 
U„ as described above for the potential calculation. The 
same parameters were also used, X = 0.2, V 0 = 1, H= 1, and 
L= 1. In addition we have let p = 1. (This just sets the overall 
scale and we are interested in relative behavior.) We evaluate 
Eqs. (56) and (57) with the coefficients T m and U n given by 
Eqs. (52) and (46). In Fig. 5 we show the vector field for the 
current density on a uniform grid of points. The direction of 
the current density, J, at a point is determined by the x and y 
components. The relative magnitude of the current density is 


y 



Fig. 5. Vector field of the current density on a uniform grid. The arrows 
indicate relative magnitude and direction. The three solid lines indicate flow 
lines and correspond to the mean path of a charge. (The flow line is every- 
where tangent to the local J.) Calculated as described in the text with \ 
= 0.2, V 0 = 1, H= 1, L= 1. 


indicated by the length. Along the top edge the magnitude of 
the current varies with position x; however, the current en- 
ters the sample with the same angle everywhere along the 
top and leaves with the same angle everywhere along the 
bottom. Notice that this implies that more charge enters the 
upper left than the upper right. Correspondingly, more 
charge leaves the lower right than the lower left. Note that 
there is an overall flow of charge from left to right. Also 
shown are three “streamlines.” These flow lines are every- 
where tangent to the local J and correspond to paths of in- 
dividual charges. Charge entering along the left edge contin- 
ues straight down to the bottom. Charge entering in the 
central portion of the top edge follows a curved path with a 
net transport in the x direction. The amount of net left to 
right deflection is the greatest for charges entering at the 
center of the top edge. 

VII. CONCLUSION 

Most textbooks of electromagnetic theory warn the reader 
that if one abandons the standard classes of boundary condi- 
tions (Dirichlet, Neumann, or “mixed"), the solubility of 
Laplace’s equation is no longer a foregone conclusion. How- 
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ever, physically plausible examples of nonstandard problems 
appear to be rather few in number. In this paper we have 
presented a solution to a nonstandard boundary value prob- 
lem that arises naturally in the context of a surprisingly 
simple Hall effect situation. 
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APPENDIX I: REDUCTION OF THE PROBLEM TO 
AN INFINITE NUMBER OF PROBLEMS OF 
STANDARD TYPE 


Because our boundary-value problem is not of standard 
type, it is not immediately obvious from a mathematical 
point of view why it should be soluble. It is therefore worth 
investigating the relation between our problem and problems 
of the standard type. To do so, we shall expand the potential 
in powers of the parameter X and then successively solve it 
to each order. Our problem then reduces to an infinite set of 
“standard" problems where the potential or its normal de- 
rivative is known on the boundaries. Essentially one solves 
the problem to some order in X and uses that solution as 
input for the next higher order. 

Consider Laplace’s equation V 2 E = 0 with the solution 
written 


V(x,y) = y ( 0 ) (x,y) + Xy ( 1 ^(x,y) + \ 2 V < 2 ) (x,y) + • • 


(62) 


where superscripts indicate the corresponding order of X for 
each term. Now the boundary conditions are V(x,±H) = 
±y 0 and dVI dx = XdV/ dy atx = 0, x = L. It is instructive to 
see how this works for the first few terms. 


±V 0 = V (0 \x,±H) + \V {1 \x,±H) 
+ \ 2 V { 2 \x,±H) + ••• 

and 

<?y (0) <?y (1) „ dv m 


(63) 


dx 


- + X 


= X 


dx 
dV (0) 


- + \~ 


dx 


x = 0,L 


dy 


- + X 


<?y (1) . dv (2) 

— + X 2 — ; + •• 


dy 


dy 


(64) 


x = 0,L 


From Eq. (63) we equate the coefficients of successive 
powers of X on each side to obtain 


y ( 0 ) (x,±//) = ±y 0 . 


V (l \x,±H) = 0, (65) 

V ( 2 ) (x,±H) = 0,... . 

Matching powers of X in Eq. (64) we see 


o , y (0) 

dx 

dV (l) 

dx 

dV (2) 

dx 


= 0, 


; = 0,L 


x = 0,L 


x = 0,L 


dV {0) 

dy 

<?y (1) 

dy 


(66) 


x = 0,L 


(67) 


1 x = 0,L 

Notice that, for each order in X, the derivative at the left and 
right boundaries, dVI dx, is now a known function of y, 
given in terms of the derivative of the solution at the previ- 
ous order. 

y" l) (x,y ) solves Laplace’s equation subject to the bound- 
ary conditions 

V (0 \x,±H)=±V 0 (Dirichlet on top and bottom), 
dV {0) 

=0 (Neumann on left and right sides). 

dx nr 

This can be done by inspection and yields V (0 \x,y) 
= Vqv/H. Note this is just the X = 0 part of our solution, Eq. 
(36).' 

y ( 1 ) (x,y) solves Laplace’s equation subject to the bound- 
ary conditions 

y (n (x,±//) = 0 (Dirichlet on top and bottom) 


dV {l) 


dx 


dv {0) 


x = 0,L 


dy 


(Neumann on left and right sides) 


1 x = 0,L 

= known function of y 
= V 0 /H. 

Solving this problem by standard methods yields 


V ,n (x,y)= S A^sinhl 

m ' = 1,3,... 


m 77 

X cos y 


(68) 


m ' 7 t 

1 L )1 

2H 



2 H 


where 


A (i)_ ~ 8V o (~ l) (m ' + 1)/2 

m' 2 v 2 cosh(m' ttL/4H) 


(69) 


(70) 


This may be compared with the order-X terms of the solution 
given above, in Eqs. (36) and (42). Denote the order-X part 
of this solution XV / |(x,y). Then from Eq. (36) we have 

Vn 


V 1 (x,y)=-\x-~ 


+ 2 

m = 1,3,. 


rO) 


x cosh 


(71) 


where the T n ] 1 are as given by Eq. (42): 


y ( 1 ) = _ 


4y n L 


7 T~Hm 2 cosh(»; ttH/L) ' 


As Eqs. (69) and (71) do not look much alike, it is worth- 
while to show that they are exactly equivalent. 
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To see this, let us perform the following expansions. Ex- 
pand sinh[w'7T(x— LI2)I2H~\ from Eq. (69) in terms of the 
functions cos (rmrx/L), on the interval where m 

= 1,3,5,... : 


sinh 



8m' LH 

IT 


2 


cosh(w' ttLIAH) lirnr \ 
m' 2 L 2 + 4m 2 H 2 C ° S | ~ X J' 


(72) 

Similarly, expanding x — L/2 from Eq. (71) on the interval 
0=£x=SL in terms of the same cos(nnrx/L) yields 

x 




Finally, expand cosh(nnry/L) from Eq. (71) on the interval 
— H^y^H in terms of the cos(m' ny/2H), where m' 
= 1,3,5,... : 


cosh — — - y =cosh 


tmrH 


1 + 


16 m H 


2 t/2 


7 T 

< i \(m r + l)/2 


X 


(- 1) 1 


, ,, m'{m' z L z + 4m z H z ) 
m = 1,3,— v 7 


m 7 t 

Xc ° s| -2H y 


(74) 


If we substitute Eqs. (73) and (74) into Eq. (71) we obtain 
64V n LH ^ ^ + m 


V\{x,y)= 


2 2 


ir 3 ~ 3 ... m'(m' Z L Z + 4m z H z ) 


Xcos 


X cos 


2 H 


(75) 



Fig. 6. Specially chosen parallelogram geometry for the Hall effect as de- 
scribed in Appendix II. The force triangle corresponding to Eq. (5) is shown 
on all four boundaries. 


Note that now the triangle of forces takes the same orien- 
tation on all four boundaries. This means that the boundary 
condition for edges 3 and 4 is now E x = 0. We must solve 
Laplace’s equation subject to the boundary conditions 


V(x,H) = V 0 (edge 1) 

(76) 

V(x, — H)= — V 0 (edge 2) 

(77) 


And if we substitute Eq. (72) into Eq. (69), we obtain an 
expression for V (1, (x,y) which is identical to the right side 
of Eq. (75). Thus Eq. (69) is indeed equivalent to the first- 
order part of Eq. (36). 

The sequence of problems beginning with Eqs. (67) and 
(68) may be continued to whatever order is necessary. This 
reduces our original problem (posed in Sec. Ill) to an infinite 
set of standard problems with “mixed" boundary 
conditions. 10 From a mathematical point of view, this ex- 
plains why our nonstandard problem should indeed have a 
stable, unique solution. From a practical point of view, how- 
ever, the solution of this problem by reduction to a sequence 
of standard problems is infinitely more tedious! 


dV 

— — =0 (edges 3 and 4). (78) 

ox 

This is not a problem of the standard type, but in this case 
the solution is trivial: 


V 0 

V(x,y) = — y 


(79) 


everywhere inside the parallelogram. The equipotentials are 
horizontal lines. 

From Eqs. (8) and (9), 

\V 0 


Jr = 


Hp( l+X 2 ) 


(80) 


APPENDIX II: ACHIEVING A UNIFORM CURRENT 
DENSITY 

The complexity of the potential problem is drastically re- 
duced if we adopt a different geometry. Let the conducting 
plate be a parallelogram, as in Fig. 6. Two edges (1 and 2) 
are parallel to the x axis. The other two edges (3 and 4) make 
an angle with the y axis that is equal to the Hall angle 0 H 
= tan _1 \=tan UB^pnq). Thus the geometry of the sample 
must be preselected to match the material properties (p and 
n) and the externally imposed magnetic field (B 0 ). 


/= ~ y ° 
y Hp{ l+X 2 ) ' 


(81) 


Thus J is everywhere constant and parallel to edge 3. The 
magnitude of J is 


V 0 

J = . (82) 

Hpsj l+\ 2 

The combination 2H^l + \ z is, of course, just the slant 
height of the parallelogram, i.e., the length of edge 3, which 
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is the actual distance that a charge must travel through the 
material. 

The total current I is obtained by multiplying J by the 
cross-sectional area of the sample (taken perpendicular to J), 
tL cos 6 H =tL!\jl + \ 2 . Thus 


tLV o 

Hp( l+\ 2 ) ’ 


(83) 


which may be compared to Eq. (61). Although the surface 
areas of the plates shown in Figs. 1 and 6 are the same 
(2/1 L), the direction of current flow in Fig. 6 produces an 
effective increase in length, and a decrease in width, of the 
conductor. 

So, another way to express the limitations of the tradi- 
tional elementary treatment is to say that it blurs the distinc- 
tion between the geometries of Figs. 1 and 6. In the limit of 
long, thin samples and low magnetic fields the two are ap- 
proximately the same. It is possible that the geometry of Fig. 
6 may have some utility in experimental applications. 
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